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Abstract 



Within the framework of the conventional QCD sum rules, we study the pion 
two-point correlation function, i f d xe iq ' x (0\TJN(x)JN(0)\^'(p)), beyond the 
soft-pion limit. We construct sum rules from the three distinct Dirac struc- 
tures, «75 ft, 175, j50-^ u q^p u and study the reliability of each sum rule. The 
sum rule from the third structure is found to be insensitive to the continuum 
threshold, S n , and contains relatively small contribution from the undeter- 
mined single pole which we denote as b. The sum rule from the 175 structure 
is very different even though it contains similar contributions from S n and b 
as the ones coming from the "fso-^q^p 1 ' structure. On the other hand, the 
sum rule from the ^75 ft structure has strong dependence on both S n and b, 
which is clearly in constrast with the sum rule for „v ( f"P v ■ We identify the 
source of the sensitivity for each of the sum rules by making specific models 
for higher resonance contributions and discuss the implication. 
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I. INTRODUCTION 



Since first introduced by Shifman, Vainshtein and Zakharov [|l|], QCD sum rule has been 
widely used to study the properties of the hadrons 0. QCD sum rule is a framework which 
connects a physical parameter to the parameters of QCD. In this framework, a correlation 
function is introduced in terms of interpolating fields, which are constructed from quark 
and gluon fields. Then, the correlation function, on the one hand, is calculated by Wilson's 
operator product expansion (OPE) and, on the other hand, its phenomenological "ansatz" 
is constructed. A physical quantity of interest is extracted by matching the two descriptions 
in the deep Euclidean region (q 2 = — oo) via the dispersion relation. The extracted value 
therefore should be independent of the possible ansatz in order to be phyically meaningful. 

The two-point correlation function with pion, 

TL(q,p) = i J d A xe^(0\T{J N (x)J N (Q)Mp)) , W 

is often used to calculate the pion-nucleon coupling, g nN , in QCD sum rules Reinders, 
Rubinstein and Yazaki calculated g n N by retaining only the first nonperturbative term in 
the OPE. Later Shiomi and Hatsuda (SH) [fj] improved the calculation by including higher 
order terms in the OPE. SH considered Eq. ([I]) and evaluated the OPE in the soft-pion limit 

(p„-°)- 

More recently, Birse and Krippa (BK) || pointed out that the use of the soft-pion limit 
does not constitute an independent sum rule from the nucleon sum rule because in the limit 
the correlation function is just a chiral rotation of the nucleon correlation function, 

11(g) = % I d*xe**(0\T[J N {x)J N (0)]\0) . (2) 

Therefore, BK considered the sum rule beyond the soft-pion limit. However, as we will 
discuss below, there seems to be mistakes in their calculation which can invalidate their 
conclusions. Thus, it is important to re-do their calculation. 

In a recent letter ||, we have pointed out that the previous calculations of the pion- 
nucleon coupling using Eq. (|1|) have dependence on how one models the phenomenological 
side; either using the pseudoscalar (PS) or the pseudovector (PV) coupling scheme. The 
two coupling schemes are equivalent when the participating nucleons are on-shell but they 
are not usually when the nucleons are off-shell. Since, in QCD sum rules, on-shell properties 
of a particle are extracted from the far off-shell point, the extracted therefore could be 
coupling-scheme dependent. Going beyond the soft-pion limit is found to be also natural in 
obtaining g n ^ independent of the PS and PV coupling schemes. In fact, we have proposed 
that, beyond the soft-pion limit, there are three distinct Dirac structures, (1) 275 ft, (2) 
^75 > (3) l^iiv(f"P v 1 eacn °f which can in principle be used to calculate g n ^. The third 
structure was found to have the common double pole structure in the phenomenological side, 
independent of the PS and PV coupling schemes. By studying this structure, we obtained the 
coupling close to its empirical value and relatively stable against the uncertainties from QCD 
parameters. Then we ask, can we get similar stable results from the sum rules constructed 
from the other Dirac structures ? If not, what are the reasons for the differences ? In 
this work, we will try to answer these questions by studying these three sum rules and 
investigating the reliability of each sum rule. 
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QCD sum rules could depend on a specific Dirac structure considered. This aspect was 
suggested by Jin and Tang in their study of baryon sum rules. They found that the chiral 
odd sum rule is more reliable due to the partial cancellation of the positive and negative- 
parity excited baryons in the continuum. Similarly here we note that the structure (1) has 
different chirality from the other two. Therefore it will be interesting to look into these sum 
rules more closely and see if similar cancellation occurs for certain sum rules. 

The paper is organized as follows. In Section [Q|, we construct three sum rules from 
the three different Dirac structures. The spectral density for the phenomenological side is 
constructed from the double pole, the unknown single pole and the continuum modeled by 
a step function. We motivate this phenomenological spectral density in Section [TTI| by using 
some effective Lagrangians for the transitions, N — > N* and N* — > N*. In Section |IV|, we 
analyze each sum rule and try to understand the differences from the formalism constructed 
in Section |TJ. A summary is given in Section M. 



II. QCD SUM RULES FOR THE TWO-POINT CORRELATION FUNCTION 

In this section, we formulate three different sum rules for the two-point correlation func- 
tion with pion beyond the soft-pion limit. For technical simplicity, we consider the correlation 
function with charged pion, 

n(g,p) = i I d A xe^ x (0\T[J p (x)J n (0)]\n + (p)) . (3) 

Here J p is the proton interpolating field suggested by Ioffe |§, 

Jp = tabciulCjftUb] 757^4 (4) 

and the neutron interpolating field J n is obtained by replacing (u,d) — > (d, u). In the 
OPE, we only keep the diquark component of the pion wave function and use the vacuum 
saturation hypothesis to factor out higher dimensional operators in terms of the pion wave 
function and the vacuum expectation value. 

The calculation of the correlator, Eq. (|J), in the coordinate space contains the following 
diquark component of the pion wave function, 

D£ = (0|<(x)4(0)|7r + (p)) . (5) 

Here, a and f3 are Dirac indices, a and a' are color indices. The other quarks are contracted 
to form quark propagators. This diquark component can be written in terms of three Dirac 
structures, 

D a J = ^(7^ 5 ) a/3 (0|rf(0)7,7 5 n(x)| 7 r + (p)) + ^(i 75 )^<0|d(0)i 7 5^)k + (p)> 

'V (75Or ^)^ ( 0|J(0) 75 ^ M (x)|7r + (p)) . (6) 



24 

Each matrix element associated with each Dirac structure can be written in terms of pion 
wave function whose first few moments are relatively well known |§. We will come back to 
the second matrix element later. For the other two elements, we need only the normalization 
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of the pion wave functions since we are doing the calculation up to the first order in p^. In 
fact, to leading order in the pion momentum, the first and third matrix elements are given 
as 0, 

(0|J(0)7 M 7 5 w(x)|7r + G»)) = iV^UPn + twist 4 term > ( 7 ) 

f tn 2 

(0|J(0)7 5 (T^M(x)|7r + (p)) = i\Pi{p ll x v -p v x^ — — — - — -. (8) 

6{m u + m d ) 

Here we have suppressed terms higher order in pion momentum. The factor v^2 is just an 
isospin factor. The twist 4 term in Eq. (]?]) comes from the second derivative term in the 
short distant expansion of the LHS. Note that in Eq. (|8]) the factor f v m\/{m u + m d) can 
be written as — (qq)/fn by making use of Gell-Mann— Oakes— Renner relation. Although 
the operator looks gauge dependent, it is understood that the fixed point gauge is used 
throughout and the final result is gauge independent. It is then interesting to note that the 
LHS of Eq. (H) can also be expanded in x such that the matrix element that contributes is 
effectively one with higher dimension, 

f 2 

(0\d(0)^ u D a u(0)\n + (p)) = iV^ip^-p^) hm " . (9) 

o(m u + m d ) 



It is now straightforward to calculate the OPE. For the ijs p 1 structure, we obtain 



q 2 \n(-q 2 ) . £ 2 hi(-g 2 ) | {^Q 2 ) | 2(qq) 2 



2tt 2 2tt 2 ' 12g 2 ' 9/ 2 g 2 



(10) 



The first three terms are obtained by taking the second term in Eq. ([])[], while the fourth 
term is obtained by taking the third term in Eq.(||) and replacing one quark propagator 
with the quark condensate. The fourth term was not taken into account in the sum rule 
studied by BK [|J but its magnitude is about 4 times larger than the third term. So there is 
no reason to neglect the fourth term while keeping the third term. The second term comes 



from the twist-4 element of pion wave function. According to Novikov et. al [ITTJ, 8 ~ 0.2 
GeV 2 . 

The phenomenological side for the 275 ft structure obtained by using the pseudoscalar 
Lagrangian takes the form, 



V2g^ N X 



N 



ni 



+ ••-. (11) 



(q 2 — m? + it) \(q — p) 2 — m? + ie] 
The dots include contributions from the continuum as well as from the unknown single 



pole terms . The latter consists of the single pole coming from iV — > iV* transition ||T0 
When the pseudovector Lagrangian is used, there is an additional single pole coming from 
iV — > iV transition ||. These single poles are not suppressed by the Borel transformation. 



1 Note that the second term in Eq. ([n]) has slightly different coefficient from BK [Q]. Ref. Q has 
the factor 5/9 instead of our factor 1/2. The difference however is small. 
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Therefore, interpretation of the unknown single pole and possibly the continuum contain 
some ambiguity due to the coupling scheme adopted. 

In principle, g n ^ has p 2 dependence as it contains the pion form factor. As one pion 
momentum is taken out by the Dirac structure, 275 ft, we take p^ = in the rest of the 
correlator as we did in the OPE side. Then the p 2 dependence of g^N can be neglected. 
Furthermore, after taking out the factor, 275 p", the rest of the correlator is a function of one 
variable, q 2 , and therefore the single dispersion relation in q 2 can be invoked in constructing 
the sum rule. Anyhow, the spectral density can be written as 

(f hen {s) = -V2g nN \ 2 N m^-5{s - m 2 ) + A' 5{s - m 2 ) + p ope {s) 9(s - S n ) . (12) 

Here the second term comes from the single pole terms whose coefficients are not known. The 
continuum contribution is parameterized by a step function which starts from the threshold, 
S n . The coefficient of the step function, p ope (s), is determined by the duality of QCD. This 
is basically the imaginary part of Eq. (0) but, because of the continuum threshold, only 
the first two terms in Eq. ( |T0| ) contribute to the coefficient. 

The parameterization of the continuum with a step function is usually adopted in the 
baryon mass sum rules. This is because each higher resonance has a single pole structure 
with a finite width. Spectral density obtained by adding up all those single poles can 
be effectively represented by a step function starting from a threshold. But in our case 
of the correlation function with pion, this parameterization for the continuum could be 
questionable. Therefore, it will be useful to construct the spectral density explicitly for 
higher resonances by employing some effective models for N* and see if the parameterization 
does make sense. This will be done in the next section. This will eventually help us to 
understand how each sum rule based on a different Dirac structure leads to different results. 

To construct QCD sum rule for the 275 ft structure, we integrate p ope (s) and p phen (s) 
with the Borel weighting factor e~ s l M and match both sides. More specifically, the sum 
rule equation after the Borel transformation is given by 

/ dse~ s/M \p ope {s) - ff hen (s)\ = . (13) 
Jo 



Using p ope (s) obtained from Eq. ( |T0D and p phen (s) in Eq. (|T2"D, we obtain 
97TN X 2 N (1 + AM 2 ) 



fn 



Jj^ e m 2 /M 2 

m 



ElM m 6 + + m 2 ( — (^g 2 ) + my ' 



2tt 2 2tt 2 V 12 \ tt / 9/2 



(14) 



Here A denotes the unknown single pole contribution, which should be determined by the 
best fitting method. Also x n = S n /M 2 and E n (x) = 1 — (1 + £ + ••• + x n /n\) e~ x . This 
expression is crucially different from the corresponding expression in Ref. 0] where the first, 
second and third terms contain the factors, E2(x n ), E\(x n ) and E (x n ) respectively. Even 
though we do not understand how such factors can be obtained, we nevertheless reproduce 
their figure by using their formula in Ref. Q and it is shown in Fig. p] (a)0. But if Eq. ( |T4|) is 



2 In plotting Figs. [T], we did not include the last term involving (qq) 2 in Eq. (|i~4|) as this term is 
new in our calculation. 
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used instead, we get Fig. [3] (b) using the same parameter set used in Ref . [§J . The variation 
scale of g n N in this figure is clearly different from the one in Fig. |l| (a). Note that some of 
their parameters are quite different from ours used in our analysis later part of this work. 
For example, 5 2 = 0.35 GeV 2 is used in Ref. 0], which is quite larger than our value of 0.2 
GeV 2 . 

QCD sum rule for the ^o iliV q ll p v structure can be constructed similarly. We have con- 
structed the sum rule for this structure in Ref. so here we simply write down the resulting 
expression, 



g wN X 2 N (l + BM 2 ) 



(qq) 
17' 



t m 2 /M 2 



M 4 £ (x w ) , 4 



12tt 2 



3 Jn \ 7T / 216 
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(15) 



Here B denotes the contribution from the unknown single pole term. Note that, since 
one power of the pion momentum is taken out by the factor, r y^a liv q^p u , we take the limit 
Pn = in the rest of the correlator as we did in the 275 p" case. In obtaining the first and 
third terms in RHS, we have used Eq. (§) while the second is obtained by taking the first 
term in Eq. (|(J) for the matrix element D^, and replacing one propagator with the quark 
condensate. The fourth term is also obtained by taking the first term in Eq. (|6|) but in this 
case other quarks are used to form the dimension five mixed condensate, (qg s o' • Qq), which 
is usually parameterized in terms of the quark condensate, m^qq). We take Wq ~ 0.8 GeV 2 
as obtained from QCD sum rule calculation ]12| . 

Now we construct QCD sum rule for the ij 5 structure. Constructing it beyond the 
soft-pion limit is more complicated as the correlator in phenomenological side has definite 
dependence on the coupling schemes. To see this, we expand the correlator for this structure 
in Pn and write 



n (g 2 ) + p ■ qli^q 2 ) + p 2 U 2 (q 2 ) + 



(16) 



Since p^ is an external momentum, the correlation function at each order of p^ can be used to 
construct an independent sum rule. Within the PS coupling scheme, the phenomenological 
correlator up to p 2 order is 



V2g nN X 



N 



p-q 



+ 



q 2 — m 2 (q 2 — m 2 ) 2 (q 2 — m 2 ) 2 



V2X 2 N p 2 dg w 



N 



q 2 — m 2 dp 2 



{p 2 = 0) 



(17) 



The dots here represent not only the contribution from higher resonances but also terms 
higher than p 2 . The last term is related to the slope of the pion form factor at p 2 = 0. 
Even though this can be absorbed into the unknown single pole term such as A or B above, 
we specify it here since this possibility is new. This correlator can be compared with the 
corresponding expression in the PV coupling scheme, 



V2g nN X 



p 2 /2 



m 



2\2 



+ 



(18) 



Note here that there are no terms corresponding to ITo and ELl. No such terms can be 
constructed from iV — > N* or N* — > N* transitions within the PV scheme. 

The single pole in Eq.flTTD survives in the soft-pion limit, which has been used by SH 
for their sum rule calculation of g n N. However, if the phenomenological correlator in PV 
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scheme is used, such sum rule cannot be constructed. Thus, going beyond the soft-pion limit 
seems to be natural for the independent determination of the coupling. However, similarly 
for n case, a sum rule can not be constructed for III. For n 2 , a sum rule can be constructed 
either in the PS or PV coupling scheme, but the residue of the double pole in Eq. ([Tj|) is 
a factor of two smaller than the corresponding term in Eq. fllTD- So the coupling-scheme 
independence can not be achieved in any of these sum rules. This is true for even higher 
orders of p M . 

A sum rule, independent of the coupling schemes, can be constructed by imposing the 
kinematical condition, 

p 2 = 2p ■ q . (19) 



With this condition, the two double pole terms in Eq. ( |TTD can be combined to yield the same 
expression as in Eq. (|18|), thus providing a sum rule independent of the coupling schemes. 
This condition comes from the on-shell conditions for the participating nucleons, q 2 = m 2 
and (q — p) 2 = m 2 , at which the physical ttNN coupling should be defined. 

The sum rule constructed with the kinematical condition, Eq. (|T9"|), is equivalent to 
consider Il 1 (g 2 )/2 + n 2 (g 2 ). This sum rule seems fine in the PS coupling scheme as there are 
nonzero terms corresponding to 111 and II2. In the OPE, the diquark component contributing 
to 275 structure is the second element of Eq. (^) which can be written in terms of twist-3 
pion wave function as |§ 

(0\d(0)n 5 u(x)\7r + (p)) = V U - / due-^ x <p p (u) . (20) 

The terms linear and quadratic in p^ in the RHS constitute the OPE correlator for II 1 and 
n 2 . Therefore, within the PS scheme, IIi and n 2 are well defined in both sides. 

Situation becomes subtle when the PV coupling scheme is employed. Before the condition 
of Eq. (|T9| ) is imposed, a sum rule can be constructed only for Il 2 as there is no Hi part in the 
phenomenological part. But after the condition, the phenomenological side has only H^ en 
which should be matched with n° pe /2 + H^ e . This seems a little awkward. Nevertheless, 
to achieve the independence of the coupling schemes, we construct a QCD sum rule for 
275 within the kinematical condition, Eq. ([19]). To be consistent with the expansion in the 
phenomenological side, we take the terms up to the order p 2 in the expansion of Eq. (P0|). 
Using the parameterization for <p p (u) given in Ref. §J, we obtained up to p 2 , 

(0|J(0)z 7 5^x)|vr + (p))= ^f 11 (l-zl-p-x-^ip-x) 2 ) . (21) 

Tfl u + Tfld \ L L J 

A different parameterization given in Ref. changes the numerical factors very slightly . 

Using the diquark component of Eq. (0), the OPE side for Hi/ 2 + n 2 is calculated 
straightforwardly. By matching with its phenomenological counterpart and taking the Borel 
transformation, we get 



(QQ) „m 2 /M 2 



7T 2 24 \ 7T 



(22) 



Here C again denotes the unknown single pole term which is not suppressed by the Borel 
transformation. This sum rule is different from the other two sum rules as its first term in 
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the OPE is negative. Each term contains very small numerical factors due to the cancellation 
between Il° pe and IlJf 6 . 

Up to now, we have presented three different sum rules from Eq. (|1]). All these sum 
rules, in principle, can be used to determine the pion-nucleon coupling constant, g n ^. We 
will discuss the reliability of each sum rule below. An alternative approach is to consider the 
nucleon correlation function in an external axial field as done in Ref ||. The nucleon axial 
charge, g^, calculated in Ref. 0, agree well with experiment. Subsequently, by using the 
Goldberger-Treiman relation, g n ^ can be also well determined. In the approach by Ref. 0, 
sum rules for — 1 is obtained by replacing some part of the OPE with the nucleon mass 
sum rule. The connection between the sum rules using Eq. (|1]) and the ones in Ref. || is 
not clear at this moment. An important observation made in Ref. || is to note that some 
(dominant) terms of the OPE correspond to a sum rule with g& = 1. This observation allows 
the construction of the sum rules for g& — 1. In our sum rules, this kind of observation is 
not possible. Also the OPE expression from Eq. (JI]) is not simply related to the sum rules 
with the external field. Therefore, Eq. (P seems to provide independent sum rules from 
the ones in Ref. ||. In future, however, further study is necessary to clarify the connection 
between these two sets of independent sum rules as it might provide important aspects for 
the nonperturbative nature of hadrons. 

III. CONSTRUCTION OF THE UNKNOWN SINGLE POLE AND THE 

CONTINUUM 

In this section, we construct the unknown single pole term and the continuum by using 
effective models for the higher resonances. This will provide a better understanding of the 
parameterization for the continuum in Eq. ([12]) and give further insights for the unknown 
single pole term. Later, this construction will help us to understand the differences between 
each sum rule based on a different Dirac structure. 

There are two possible sources for the unknown single pole term and the continuum. 
One is from the transition, N — > N* and the other is from the transition, N* — > N*. Of 
course, as we pointed out in Ref. ||, there could be additional single pole of nucleon coming 
from N —>■ N, which however, in the first order of the pion momentum, appears only in 
the sum rule for the ijs $ structure within the PV coupling scheme. First we avoid such 
possibility by constructing effective models within the PS coupling scheme. Later we will 
discuss the case with the PV coupling scheme. Moreover, we will discuss only the two Dirac 
structures, 275 p" and 75a iiv q i *"p v '. The correlator for the 275 structure with the kinematical 
condition of Eq. flTP| ) takes almost the same form as the one for the 75a ^q^p" structure. A 
slight difference is the appearance of terms containing the derivative of pion form factor as 
indicated in Eq. (|T7|) . Note that this difference is only specific to the PS coupling scheme. 
As the form factor is a smooth function of p 2 around p 2 = 0, this difference is not expected 
to be crucial. 

Within the PS coupling scheme, N —>■ N* contributions to the correlator, Eq. (^), can 
be constructed by using the effective Lagrangians for the positive and negative (tp-) 
parity resonances, 



s 



gnNN-^T ■ Kip- - gnNN.tp-iT • 7T^ . (23) 

The nucleon field is denoted by ip here. These terms contribute to the correlator because 
the nucleon interpolating field can couple to the positive and negative parity resonances via, 

(0\J N \N + (k,s)) = X + U(k,s); (0\J N \N_(k,s)) = \- l5 U(k,s) , (24) 

where U (k, s) denotes the baryon Dirac spinor and A± indicates the coupling strength of the 
interpolating field to each resonance with specified parity. 
The 75a iiv q ll "p v structure of the correlator takes the form, 

2AjvA_<77rArjv_ 2Aat A + (? 7r jvjv + /__s 

+ 712 ™2W„2 Z^2 \ ' \ Zb ) 



(q 2 — m 2 )(q 2 — m?) (q 2 — m 2 )(q 2 — m\) ' 

which can be compared with the ^75 $ structure 

2X N X_g vNN _ (m_ - m) _ 2X N X + g 7tNN+ (m + + m) 
(q 2 — m 2 )(q 2 — m 2 _) (q 2 — m 2 )(q 2 — m 2 + ) 

By separating as 

1 



(26) 



(q 2 — m 2 )(q 2 — m±) m± — in 2 



q 2 — m 2 q 2 — 171% 



(27) 



we can see that the transitions, iV — > A^*, involve the two single poles, one with the nucleon 
pole and the other with the resonance pole. The former constitutes the unknown single pole 
as it involves the undetermined parameters, A± and g n NN±- I n the latter, the finite width 
of the resonances can be incorporated by replacing m± — > m± — «r±/2 in the denominator. 
Then when it is combined with other such single poles from higher resonances, it produces 
the spectral density which can be parameterized by a step function as written in Eq. (|12|). 
This also implies that the continuum threshold, S n , does not need to be different from the 
one appearing in the usual nucleon sum rule. 

It is now easy to obtain the spectral density for the two Dirac structures by incorporating 
the decay width of the resonances. For the ^75 $ structure, we have 

/ \ o ( ^+gnNN+ . X-g nNN _ \ . 2 

ps{s) = 2 i H ■ X N 5{s - m ) 

y m + — m m_ + m J 

+ 2X ^ X ^nn + _ 2X N X-g vNN _ 

m + — m ' m_ + m 

where 

G(s, m±) = - ? ?,f_f 2p2 • (29) 

Note that the contribution from the positive-parity resonance is enhanced by the factor 
1 / (m+ — m) while the one from the negative-parity resonance is suppressed by the factor 
l/(m_ + m) . Similarly for the ^c^q^p" structure, we obtain 
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ps{s) = 2 — = 1 H = \ N 6{s - m ) 



mi — mr m_ — m" 



mi — rrr mi — 



Note that the superficial relative sign between the positive- and negative-parity resonances 
are opposite to that in Eq. (|28|). It means, depending on the relative sign between X + g 7T NN + 
and \-g n NN_, the two contributions add up in one case or cancel each other in the other 
case. In other words, we can say something about the coefficients of S(s — m 2 ) and G(s, m±), 
by studying the sensitivity of the sum rules to the continuum or to the single pole. 

Additional contribution to the continuum may come from N* — > N* transitions. For the 
off-diagonal transitions between two parities, N + — ► N- and i\L — > N + , we use the effective 
Lagrangians, 

gnN+N.i'+ir ■ 7T^_ - gnN+N.lp-iT • TT1p + , (31) 

to construct the correlator. These off-diagonal transitions lead to the spectral density of 

Pod{s) oc A_A+[G(s,m + ) - G(s,m_)] , (32) 

which is therefore suppressed by the cancellation between the two parity resonances. 

For the diagonal transitions, iV + — > iV + and iV_ — ► A^_, we use the effective Lagrangians, 

g-nN+N+^+ilsT ■ TT1p + ; g-xN-N.^-ilbT ■ 7T?/>- . (33) 

These diagonal transitions produce only the double pole for the correlator, 1/ (q 2 — m 2 ± + 
im±r±) 2 , which is then translated into the spectral density, 

-m±g*N ± N ± >?±j-G{s,m±) for ry 5 jf 
Pd(s) ~ { ^ (34) 

±g7rN ± N ± >?±£G(s,m±) for ^a^p v . 



First note that, because of the derivative, each spectral density has a node at s = m|, posi- 
tive below the resonance and negative above the resonance. Then under the integration over 
s, the spectral density from the double pole is partially canceled, leaving attenuated contri- 
bution coming from the s dependent Borel weight. Indeed, one can numerically check that, 
for the Roper resonance, Jdse- s / M2 G(s , m + ) is always larger than / dse S//m2 dG(s,m + )/ds 
for M 2 > 0.7 GeV 2 and the cancellation is more effective as M 2 increases. In general, the 
continuum contributes more to a sum rule for larger M 2 . Hence the double pole is more 
suppressed than the single pole in the region where the continuum is large. Further sup- 
pression of the double pole continuum can be observed, for example, by comparing the first 
equation of Eq. ( |34"D with Eq. (p8|). Even if one assumes that gnNN + ^N^+ ~ fiVTV+iv+A^ , 



3 The nucleon interpolating field, Jn, is constructed such that it couples strongly to the nucleon 
but weakly to excited states. Therefore, A+ is expected to be smaller than Aat. This assumption, 
therefore, may be regarded as assuming strong coupling to the excited baryon. 
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then Eq. fl28|) has the enhancing factor of 1/ (m + — m) while the first equation in Eq. fl34|) 
contains only m + . Thus, the double pole contribution is much suppressed than the single 
pole, which can be checked also from numerical calculations. The similar suppression can 
be expected for the second equation in Eq. ([34]). Therefore, we expect that the continuum 
mainly comes from the single pole of l/(q 2 — m\ + im±Y±) which is generated only from the 
N — > N* transitions. This will justify the "step-like" parameterization of the continuum as 
given in Eq. flT2"|). 

Now we discuss the case with the PV coupling scheme. We use the following Lagrangians 

9nN+N+ 7 ™ , g-TrN-N- j o/t , 

-^M +7W - • + (H.C.) , 

m + m + 



m_ — m 

g-irN+N- 

m_ — m + 



^_ 7 57^ • + (H.C.) . (35) 



These effective Lagrangians in the PV scheme are constructed such that the action is the 
same as the PS case when the resonances are on-shell . In this case, complications arise 
from the possible single pole term coming from iV — > N |J contribution which was absent 
in the PS scheme. This also means that there could be additional single poles coming from 
N — > iV* and N* — > iV* transitions. Note, this kind of complication arises only in the 275 
case. That is, for the "fscr^q^p 1 ' case, we have the same spectral density as given in Eq. (|30|). 

As we mentioned above, because the double pole type contribution, l/(g 2 — m±+im±r±) 2 , 
to the continuum is suppressed, only single poles are important in constructing the spectral 
density for the unknown single pole and the "step-like" continuum. To construct the single 
poles, we consider all possibilities, N — > N, N — > iV* and iV* — > iV*. The coefficient of 
\n8(s — m 2 ), namely the unknown single pole term for the 275 structure, can be collected 
from A" — > N and A" — > N* transitions, 



g-irN^N _ 2 m ( ^+g-irNN+ ^~girNN. 

2m \ ml — m 2 ml — m 2 



- 2m [—2 Z^ + —2 =3 • ( 36 ) 



Compared with the corresponding term in Eq. fl3T|), this term is differed by the first term 
associated with A" — > N. The second and third terms are the same except for the overall 
factor, —2m. 

Also the continuum contributions are collected from the terms containing l/(q 2 — m^) 
in the correlator. We thus obtain the spectral density for the continuum, 

2m+AArA + 2m + A + A 




+ gnNN+ 9 9 gnN+N 



m + — m z m_ — 



, ,._ 2m_AivA_ 2m_A+A_ . 

+ [g-TrN-N-- g-n-NN. 9 9 g-rrN+N- 9 T ^\ s i m -) • [< i <) 



9 9 y-7riv + iv_ 9 

m_ — m mi — m 
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IV. RELIABILITY OF QCD SUM RULES AND POSSIBLE INTERPRETATION 



In section O, we have constructed three sum rules, each for the ft, the and the 
j5<J fJiU q fl p lJ structures beyond the soft-pion limit. Ideally, all three sum rules should yield the 
same result for g n ^. In reality, each sum rule could have uncertainties due to the truncation 
in the OPE side or large contributions from the continuum. Therefore, depending on Dirac 
structures, there could be large or small uncertainties in the determination of the physical 
parameter. This can be checked by looking into the Borel curves and seeing whether or 
not they are stable functions of the Borel mass. In the QCD sum rules for baryon masses, 
the ratio of two different sum rules is usually taken in extracting a physical mass without 
explicitly checking the stability of each sum rule. As pointed out by Jin and Tang j7|, this 
could be dangerous. In this section, we will demonstrate this issue further by considering 
three sum rules provided in section ||. 

In Eqs. (|nj), and flT5p, LHS can be written in the form, c + bM 2 . The parameter c 
denotes the same quantity, i.e. gvjvAjv, but b could be different in each sum rule. We can 
determine c and b by fitting RHS by a straight line within the appropriately chosen Borel 
window. Usually, the maximum Borel mass is determined by restricting the continuum 
contribution to be less than, say, 30 ~ 40 % of the first term of the OPE and the minimum 
Borel mass is chosen by restricting the highest dimensional term of the OPE to be less than, 
say 10 ~ 20 % of the total OPE. These criteria lead to the Borel window centered around the 
Borel mass M 2 ~ 1 GeV 2 . Further notice that c determined in this way does not depend on 
the PS and PV coupling schemes while the interpretation of b could be scheme-dependent. 

In the analysis below, we use the following standard values for the QCD parameters, 

(qq) = -(0.23 GeV) 3 ; (^£ 2 ) = (0-33 GeV) 4 , 

5 2 = 0.2 GeV 2 ; m 2 = 0.8 GeV 2 . (38) 

Uncertainties in these parameters do not significantly change our discussion below. For the 
nucleon mass m and the pion decay constant /„-, we use their physical values, m = 0.94 GeV 
and U = 0-093 GeV. 

In Figure |2] (a), we plot the Borel curves obtained from Eqs.(|I4]), ( |2"2"D and flT5|). The thick 
solid line is from Eq.flTBD, the thick dot-dashed line from Eq.(|22"D and the thick dashed line 
from Eq.([i4|). In all three curves, we use S n = 2.07 GeV 2 corresponding to the mass squared 
of the Roper resonance. To check the sensitivity on S^, we have increased the continuum 
threshold by 0.5 GeV 2 and plotted in the same figure denoted by respective thin lines. 

In extracting some physical values, one has to fit the curves within the appropriate Borel 
window using the function c + bM 2 . The unknown single pole term, 6, is represented by the 
slope of each Borel curve. The intersection of the best fitting curve with the vertical axis 
gives the value of c. Figure |2| (b) shows the best fitting curves within the Borel window, 
0.8 < M 2 < 1.2 GeV 2 . This window is chosen following the criteria mentioned above. But 
as the Borel curves are almost linear around M 2 ~ 1 GeV 2 , the qualitative aspect of our 
results does not change significantly even if we use the slightly different window. 

The r )^cr liu q^p v sum rule yields c ~ 0.00308 GeV 6 . To determine g n N, the unknown 
parameter X N needs to be eliminated by combining with the nucleon odd sum rule [|J. 
According to the analysis in Ref. ||, this sum rule yields g T N ~ 10 relatively close to its 
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empirical value. As can be seen from the thin solid curve which is almost indistinguishable 
from the thick solid curve in Fig. ^| (b), this result is not sensitive to the continuum threshold, 
S n . Also note from Table | that the unknown single pole term represented by b is relatively 
small in this sum rule. 

The result from the 1^5 sum rule is c ~ —0.0003 GeV 6 , which is obtained from linearly 
fitting the thick dot-dashed curve in Fig. |2] (a). Even though the thin dot-dashed curve is 
almost indistinguishable from the thick dot-dashed curve, the best fitting value for c with 
S n = 2.57 GeV 2 is about 50 % smaller than the one with SV = 2.07 GeV 2 . This is because 
the total OPE strength of this sum rule is very small. The negative value of c indicates 
that g n N is negative. Also the magnitude of this is about a factor of ten smaller than 
the corresponding value from the 75a ilv <f"P v sum rule. When this result is combined with 
the nucleon odd sum rule, then the extracted ttNN coupling would be a lot smaller than 
its empirical value and therefore it can not be acceptable as a reasonable prediction. As 
we discussed in Section ||, the problem might due to the kinematical condition, Eq. flTSp. 
Though we have introduced this condition in order to achieve the independence from the 
coupling scheme employed, this condition inevitably combines two independent sum rules, 
111 and II2 in Eq. fll6|), which reduces the OPE strength. This reduction makes the 275 
sum rule less reliable because of the cancellation of the main terms. Nevertheless, this 
study shows that one could get a totally different result depending on how the sum rule is 
constructed. 

For the 275 sum rule, the Borel curve around M 2 ~ 1 GeV 2 is almost a linear function 

of M 2 . By linearly fitting the thick dashed curve (5 W = 2.07 GeV 2 ), we get c 0.00022 

GeV 6 . But with using S v = 2.57 GeV 2 , we obtain c ~ -0.0023 GeV 6 , a factor of ten 
larger in magnitude. Thus, there is a strong sensitivity on S w which changes the result 
substantially. Again c is negative in this sum rule, indicating that g n ^ is negative. The sign 
of this result however depends on the Borel window chosen. Restricting the Borel window 
to smaller Borel masses, the extracted c becomes positive though small in magnitude. The 
slope of the Borel curve is also large, indicating that there is a large contribution from the 
undetermined single pole terms. The thin dashed curve ( for S n = 2.57 GeV 2 ) in Fig. || (b) 
is steeper than the thick dashed curve (for S n = 2.07 GeV 2 ). In a sum rule, the larger 
continuum threshold usually suppresses the continuum contribution further. Since a more 
steeper curve is expected as we further suppress the continuum, this sum rule contains very 
large unknown single pole terms. This provides a very important issue which should be 
properly addressed in the construction of a sum rule. The unknown single pole terms could 
be small or large depending on a specific sum rule one considers. 

From the three results, we showed that the extracted parameter, here c, could be totally 
different depending on how we construct a sum rule. Even the sign of the parameter is not 
well fixed. Certainly the r )^ > o lliV q^p v sum rule has nice features, such as small contributions 
from the continuum and the unknown single pole. And when it is combined with the nucleon 
odd sum rule, it provides g^N reasonably close to its empirical value ||. But the other sum 
rules do not provide a reasonable or stable result. It is not clear if this is due to the lack of 
convergence in the OPE or due to the limitations in the sum rule method itself. To answer 
such questions, it would be useful to analyze the OPE side further. However our analysis 
raises an issue whether or not a sum rule based on one specific Dirac structure is reliable. 

Still, regarding the sensitivity of S w and the unknown single pole contribution, we can 
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provide a reasonable explanation based on effective model formalism developed in Section JU\ . 
Results from the two sum rules, 275 and r ysG tw q tt p v structures, share similar properties. As 
can be seen from Table |, for the 275 sum rule, the extracted c is —0.00033 GeV 6 when 
S w = 2.07 GeV 2 is used. For S n = 2.57 GeV 2 , c = -0.00016 GeV 6 . So the difference is 
0.00017 GeV 6 . This difference is close to the difference from the , j5<J llu q ll p L/ case. Further- 
more, the magnitude of b is relatively close in the two sum rules. These common behaviors 
of the two sum rules are expected because, as we briefly mentioned in section JTJ, their 
phenomenological structures for the continuum and the unknown single poles are almost the 
same except for the possible small term containing the derivative of the pion form factor. 
[See Eq. (0).] The similar slope and the similar contribution from S n are actually related 
as can be seen from Eq. (^). In Eq. (|30|), the terms corresponding to the unknown single 
poles have the same relative sign between the positive- and negative-parity resonances as 
the terms corresponding to the continuum. If we assume that the sign of \ + g n NN + is op- 
posite to that of \-.g n NN-, then there is a cancellation between the two resonances. Thus, 
with this sign assignment, we expect both terms, unknown single pole and the "step-like" 
continuum, contribute less to the sum rules. This is what Fig. ^] indicates. As Eq. ( p0|) is 
independent of the coupling schemes, this explanation is valid even for the PV case. The 
sign assignment, within the PS coupling scheme, also explains the large slope and strong 
sensitivity of S n in the ^75 ^ sum rule. From Eq. (^) with the sign assignment, negative- 
and positive-parity resonances add up for the undetermined single pole and the continuum, 
yielding large contribution to the two. 

This explanation for the 275 j6 sum rule can be changed for the PV coupling scheme. 
For the case with the undetermined single pole, as can be seen from Eq. (p6[), resonances 
with different parities cancel each other also for the 275 p" case under the sign assignment 
introduced above. However, there is an additional single pole coming from N — > N which 
could explain the large slope. Its contribution to A in Eq. (|14D can be calculated to be 
—l/2m. In terms of mag nitude, it contributes 50% of the LHS at M 2 ~ 1 GeV 2 with the 
opposite sign from the first term. Since c is negative as we showed in Table |, g^N is also 
negative. Since b ~ gy/vA the unknown single pole term is positive which can explain the 
large and positive slope in this sum rule. As for the continuum, Eq. fl3"T|) shows that there are 
other contributions associated with N* — > N* whose magnitudes can not be estimated. Even 
though we can not say that the large continuum only comes from adding up the positive- 
and negative-parity resonances, this is not contradictory to the sign assignment for g n NN + ^+ 
and g n NN-^~- Note however that the negative sign of c is not firmly established in this sum 
rule for the 275 p 1 structure as there is a possibility that c can be positive for different Borel 
window chosen. In this case, the positive and large slope of the Borel curve can not be well 
explained within the effective model. 

Nevertheless, our study in this work, though it was specific to the two-point nucleon 
correlation function with pion, raises important issues in applying QCD sum rules in cal- 
culating various physical quantities. Most QCD sum rule calculations are performed based 
on a specific Dirac structure without justifying the use of the structure. As we presented in 
this work, a sum rule result could have a strong dependence on the specific Dirac structure 
one considers. This dependence is driven by the way how the sum rule is constructed or 
by the difference in the continuum contributions or the unknown single pole terms. The 
continuum and the unknown single pole terms are large in some case while they are small 
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in other cases. 



V. SUMMARY 

In this work, we have presented three different sum rules for the two-point correlation 
function with pion, i J d A xe iq ' x {Q\T Jn{x)Jn{^)\^{j>)) i beyond the soft-pion limit. The PS 
and PV coupling scheme independence has been imposed in the construction of the sum 
rules. We have corrected an error in the previous sum rule in Ref. jp and found that the 
sum rule contains large contribution from the unknown single pole, b, and the continuum. 
On the other hand, the sum rules for ij^ and 75a '^vQ^p 11 structures share similar properties, 
relatively similar contributions from the continuum and the unknown single pole. By making 
specific models for higher resonances, we have explained how the latter two sum rules are 
different from the ^75 $ sum rule. Within the PS coupling scheme, the difference can be well 
explained by the cancellation or addition of the positive- and negative-parity resonances in 
higher mass states. Within the PV coupling scheme, the large slope of the Borel curve in 
the 275 j6 sum rule can be attributed to the single pole coming from N — > N transition 
even though this explanation is limited to the case with negative value of g n N- The value 
of c extracted from the 275 and 75a f _ lu q fJ 'p u sum rules are different. For the 275 sum rule, in 
order to eliminate the coupling scheme dependence, we need to impose the on-mass-shell 
condition before the matching the OPE and phenomenological correlators. Then a significant 
cancellation occurs and it makes the 275 sum rule less reliable. We have stressed that in the 
construction of a sum rule, a care must be taken. 
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TABLES 



TABLE I. The best-fit values for the parameters c and b obtained within the Borel window 
0.8 < M 2 < 1.2 GeV 2 . The numbers in parenthesis are obtained when S w = 2.57 GeV 2 is used. 





c (GeV 6 ) 


b (GeV 4 ) 


*75 i> 


-0.00022 (-0.0023) 


0.011 (0.0145) 




-0.00033 (-0.00016) 


-0.00183 (-0.0021) 




0.00308 (0.002906) 


0.00257 (0.0029) 
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FIGURES 



FIG. 1. The result of Birse and Krippa's sum rule for g n N is shown in Figure (a). The 
solid line is obtained by eliminating the unknown single pole term using the differential operator 
1 — M 2 d/dM 2 and the dashed line is obtained simply by neglecting the unknown single pole term. 
Figure (b) is similarly obtained but after correcting the factors in the treatment of the continuum. 

FIG. 2. (a) Borel curves obtained from the three different sum rules. The thick solid line 
(thick dashed line) is for the -y^a^ q^Pv («75 tf) structure. The thick dot-dashed line is for the 275 
structure. Corresponding thin lines are obtained when S w = 2.57 GeV 2 is used, (b) The curves 
obtained by linearly fitting the Borel curves within the range, 0.8 < M 2 < 1.2 GeV 2 . The thin solid 
lines and thin dot-dashed lines are almost indistinguishable from the corresponding thick lines. 
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